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Parity -time-symmetric  whispering-gallery 
microcavities 

Bo  Peng1*,  §ahin  Kaya  Ozdemir1*4,  Fuchuan  Lei1,2,  Faraz  Monifi  \  Mariagi  ovanna  Gianfreda3,4, 
Gui  Lu  Long2,5,  Shanhui  Fan6,  Franco  Nori7,8,  Carl  M.  Bender3  and  Lan  Yang1* 


Optical  systems  combining  balanced  loss  and  gain  provide  a  unique  platform  to  implement  classical  analogues  of  quantum 
systems  described  by  non-Hermitian  parity-time  (PT)-symmetric  Hamiltonians.  Such  systems  can  be  used  to  create  synthetic 
materials  with  properties  that  cannot  be  attained  in  materials  having  only  loss  or  only  gain.  Here  we  report  PT-  symmetry 
breaking  in  coupled  optical  resonators.  We  observed  non-reciprocity  in  the  PT-symmetry-breaking  phase  due  to  strong  field 
localization,  which  significantly  enhances  nonlinearity.  In  the  linear  regime,  light  transmission  is  reciprocal  regardless  of 
whether  the  symmetry  is  broken  or  unbroken.  We  show  that  in  one  direction  there  is  a  complete  absence  of  resonance  peaks 
whereas  in  the  other  direction  the  transmission  is  resonantly  enhanced,  a  feature  directly  associated  with  the  use  of  resonant 
structures.  Our  results  could  lead  to  a  new  generation  of  synthetic  optical  systems  enabling  on-chip  manipulation  and  control 
of  light  propagation. 


Parity-time  (PT)-symmetric  quantum  Hamiltonian  systems 
have  attracted  increasing  attention  since  it  was  shown  that 
the  eigenvalues  of  non-Hermitian  Hamiltonians  H 

can  be  entirely  real  iff  hey  respect  PT  symmetry1,  PTH  =  HPT. 
PT-symmetric  systems  are  open  physical  systems  having  balanced 
absorption  (loss)  and  amplification  (gain).  Remarkably,  such 
systems  can  exhibit  a  phase  transition1-3  (spontaneous  PT- 
symmetry  breaking)  iff  he  parameter  controlling  the  degree  of 
non-Hermiticity  exceeds  a  critical  value.  Beyond  this  threshold 
the  eigenvalues  become  complex  even  though  PTH  =  HPT  still 
holds.  PT  symmetry  has  been  studied  experimentally4-14  and 
theoretically15-25  in  various  physical  systems,  with  optics  emerging 
as  the  most  versatile  platform  to  explore  PT-symmetric  applications. 

So  far,  experiments  in  PT-symmetric  optics  have  been  limited 
to  waveguides4-8  in  which  resonances  play  no  role.  Loss-induced 
transparency4,  power  oscillations  violating  left-right  symmetry5, 
PT-synthetic  photonic  lattices7,  and  unidirectional  invisibility7,8 
have  been  demonstrated,  but  other  phenomena  such  as  non¬ 
reciprocal  light  transmission6,15,16,26,27,  coexisting  coherent  perfect 
absorption  and  lasing28-30,  photonic  analogues  of  topological 
insulators31  and  exceptional  points  in  lasers32  are  yet  to  be 
realized.  These  could  benefit  significantly  from  resonance  structures 
exhibiting  PT  symmetry.  In  particular,  optical  non-reciprocity27,33-38 
has  been  a  long- sought- after  goal  in  the  study  of  PT  symmetry 
because  ofi  ts  substantial  device  implications  for  optical  information 
processing.  Here  we  report  PT-symmetry  breaking  in  a  system 
of  two  directly  coupled  on-chip  optical  resonators,  and  provide 
direct  proof  of  non-reciprocity  in  PT-symmetric  optics.  Our  record 
low  power  of  1  pW,  to  observe  nonlinearity-induced  time-reversal 
symmetry  breaking  for  non-reciprocal  light  transmission,  could  not 
be  achieved  without  a  resonant  structure. 


Our  system  consists  off  wo  directly  coupled  microtoroidal 
whispering-gallery-mode  resonators39,40  (WGMRs),  each  coupled 
to  a  different  fibre-taper  coupler  (Fig.  la-c).  This  system  is  PT 
symmetric  because  under  parity  reflection  P  the  WGMRs  become 
interchanged  and  under  time  reversal  T  loss  becomes  gain  and  gain 
becomes  loss.  The  first  microtoroid  (pRj)  is  an  active  resonator 
made  from  Er3+ -doped  silica40-42;  the  second  microtoroid  (pR2)  is 
a  passive  (no-gain-medium)  resonator  made  from  silica  without 
dopants40.  Gain  in  pR!  was  provided  in  the  1,550  nm  wavelength 
band  by  optically  pumping  Er3+  ions  with  a  pump  laser  in  the 
1,460  nm  band.  The  Q-factors  of  pR,  and  pR2  in  the  1,550  nm  band 
were  3.3  X  106  and  3  X  107,  respectively,  and  pR,  had  a  Q-factor 
of  2.4  X  106  in  the  1,460  nm  band  (Fig.  Id).  The  microtoroids 
were  fabricated  at  the  edges  off  wo  separate  chips  placed  on 
nanopositioning  systems  to  control  precisely  the  distance  and  hence 
the  coupling  between  the  microtoroids  (Supplementary  Sections  A 
and  Bl).  We  mediated  the  coupling  between  pR,  and  pR2  in  the 
1,550  nm  band  by  controlling  the  detuning  between  their  resonant 
wavelengths  through  the  tuning  off  he  resonance  wavelength  of  pR2 
via  the  thermo-optic  effect  of  silica.  There  was  no  coupling  between 
the  resonators  in  the  1,460  nm  band;  thus,  the  pump  existed  only 
in  pR,.  Compensation  of  the  pR,  losses  in  the  1,550  nm  band  with 
the  optical  gain  provided  by  Er3+  was  confirmed  by  the  narrowing 
of  resonance  linewidth  with  increasing  pump  power  (Fig.  le)  and 
by  the  emergence  of  a  strong  resonance  peak  (Fig.  If)  due  to  the 
amplification  of  a  very  weak  probe  by  the  gain. 

We  conducted  two  sets  of  experiments  using  the  apparatus  in 
Fig.  1.  The  first  set  determined  the  broken  and  unbroken  PT  phases 
in  coupled  microcavities  with  balanced  loss  and  gain  as  a  function 
of  the  coupling  strength  K,  unlike  previous  experiments  where  K 
was  fixed  and  the  gain  and  loss  ratio  was  varied.  We  studied  the 
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Figure  1 1  PT- symmetric  whispering-gallery -mode  microcavities,  a,  A  system  consists  of  two  directly  coupled  whispering-gallery-mode  resonators 
(W  GMRs )  and  two  fibre-taper  waveguides.WG1:fibre-taper  waveguide  with  ports  1  and  2.  WG2:  fibre-taper  waveguide  with  ports  3  and  4.  pRp  active 
Er3+  -doped  silica  microtoroid.  pR2:  passive  silica  microtoroid.  Pin:  pump  laser  in  1460  nm  band  to  excite  Er3+  ions  that  provide  gain  in  1,550  nm  band.SiN: 
probe  light  (signal)  in  1,550  nm  band.  b,c,  Top  and  side  views  of  the  coupled  resonators,  d,  Transmission  spectra  showing  resonance  of  pR2  at  1,533.8  nm  (i) 
and  resonances  of  piR  i  at  1,533.8nm  (ii)  and  at  1417.0  nm  (iii).  e,  Gain  provided  by  Er3+  ions  in  pRi  leads  to  narrower  and  deeper  resonances  as  pump 
power  (gain)  is  raised,  implying  an  increasing  Q-factor.  f,  Weak  probe  light  is  amplified  when  coupled  to  pRi  together  with  the  pump  light.  Inset  shows  that 
without  the  weak  Sin  there  is  no  resonance  enhancement. 


system  using  only  the  waveguide  (WG1)  coupled  to  pR^  The  pump 
and  the  weak  probe  lasers  were  input  at  port  1  and  the  output 
transmission  spectra  were  monitored  at  port  2  in  the  1,550  nm 
band.  Without  the  pump,  the  coupled-resonator  system  acted  as 
a  passive  photonic  molecule43  characterized  by  two  supermodes 
whose  spectral  distance  increases  with  increasing  k  as  seen  in 
Fig.  2a  (k  decreases  exponentially  with  increasing  distance  between 
pRj  and  pR2;  see  Supplementary  Fig.  5).  This  system  became  PT 
symmetric  when  pR!  was  optically  pumped  to  provide  gain  and 
pR2  had  a  balanced  loss.  At  fixed  gain-loss  ratio,  we  monitored  the 
output  port  as  a  function  of  k  and  observed  the  PT  phase  transition 
at  the  threshold  coupling  strength  kpt  (Fig.  2a, b).  For  k/  kpt  <  1, 
the  system  is  in  a  broken- symmetry  phase,  as  seen  in  both  the 
coalescence  of  the  real  parts  of  the  eigenfrequencies  (Fig.  2a)  and 
the  non-zero  difference  in  their  imaginary  parts  (Fig.  2b).  As  k/  kpt 
approaches  1  from  below,  the  difference  in  the  imaginary  parts 
of  the  eigenfrequencies  decreases  and  their  real  parts  bifurcate 
(mode-splitting). 

Next,  we  chose  two  different  WGMs  with  Q-factors  2.0  x  107 
and  3.0  x  107  in  pR2  and  adjusted  the  pump  power  so  that  the  loss- 
gain  ratio  was  nearly  balanced.  We  observed  the  transition  from  the 
broken  to  unbroken  phase  occurring  at  different  coupling  strengths 
for  modes  with  different  Q,  that  is,  different  initial  loss  (Fig.  2c, d). 
The  PT  phase  transition  occurs  at  higher  kpt  for  lower  Q-factors. 

The  PT  phase  transition  can  be  understood  intuitively  as  follows. 
If  the  coupling  between  the  resonators  is  weak,  the  energy  in 
the  active  resonator  cannot  flow  fast  enough  into  the  passive 
resonator  to  compensate  the  absorption.  Thus,  the  system  cannot 
be  in  equilibrium  and  the  eigenfrequencies  are  complex,  implying 
exponential  growth  or  decay.  Flowever,  if  the  coupling  strength 
exceeds  a  critical  value,  then  the  system  can  attain  equilibrium 
because  the  energy  in  the  active  resonator  can  flow  rapidly  enough 
into  the  passive  one  to  compensate  the  dissipation. 

In  our  experiments  the  frequency  bifurcation  (splitting)  is  not 
in  orthogonal  directions  (Fig.  2)  as  would  be  expected  for  ideal 


systems  with  exactly  balanced  gain  and  loss.  Instead,  the  bifurcation 
is  smooth  and  the  degree  of  smoothness  (how  much  the  system 
deviates  from  the  exactly  balanced  case)  depends  on  the  pump 
power.  To  understand  the  origin  of  this  behaviour,  we  revisited  the 
equations  of  motion  for  coupled  oscillators,  which  showed  that  for 
unbalanced  gain  and  loss,  the  eigenfrequencies  are  never  exactly 
real44'45.  Instead,  there  is  a  region  of  k  where  the  difference  in 
imaginary  parts  is  large  but  the  difference  in  real  parts  is  small  (but 
non- zero),  and  a  second  region  where  the  difference  in  imaginary 
parts  is  small  but  non-zero,  and  the  splitting  is  large.  In  practical 
implementations  it  is  impossible  to  balance  the  loss  and  gain  exactly, 
so  the  mathematical  prediction  of  a  smooth  bifurcation  is  physically 
realistic  and  consistent  with  our  experiments  (Fig.  2). 

A  linear  static  dielectric  system,  even  with  gain  and  loss, 
cannot  have  non-reciprocal  response27'33,38.  However,  a  nonlinear 
system  can  exhibit  strong  non-reciprocity.  In  our  second  set  of 
experiments  we  tested  this  in  our  PT-symmetric  system  (Fig.  la 
with  transmission  from  input  port  1  (4)  to  output  port  4(1)  defined 
as  the  forward  Ti->4  (backward  T4^i))  and  demonstrated  strong 
non-reciprocal  light  transmission  associated  with  nonlinearity 
enhancement  induced  by  PT-symmetry  breaking. 

We  first  monitored  the  output  spectra  at  port  1  as  the  power  of 
input  probe  at  port  4  was  varied  while  the  system  was  in  broken- 
or  unbroken- symmetry  phases.  A  clear  nonlinear  response  was 
observed  in  the  symmetry- broken  phase  in  contrast  to  the  linear 
response  in  the  unbroken  phase  (Fig.  3a).  At  low  power  levels  where 
the  input-output  relation  was  linear,  the  system  was  reciprocal 
in  both  the  broken-  and  unbroken- symmetry  phases  (Fig.  3b, c). 
Thus,  we  have  direct  experimental  clarification  of  reciprocity  in  PT- 
symmetric  systems;  PT  symmetry  alone  is  not  sufficient  for  non¬ 
reciprocal  light  transmission.  As  we  increased  the  input  power,  the 
system  remained  in  the  linear  regime  for  the  unbroken- symmetry 
phase,  whereas  the  input-output  relation  became  nonlinear  in 
the  broken  phase  (Fig.  3a).  These  results  indicate  nonlinearity 
enhancement  (that  is,  lower  threshold  for  nonlinearity)  in  the 
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Figure  2  |  PT- symmetry  breaking  in  coupled  WG  M  microresonators,  a-d,  Experimentally  obtained  real  Re(co)  ( a,c)  and  imaginary  Im(cj)  (  b,d)  parts  of  the 
eigenfrequencies  of  supermodes  as  a  function  of  coupling  k/kpt  (Supplementary  Sections  B2  and  B3  and  Supplementary  Fig.  6).  (Pq  is  the  frequency  of  the 
resonances  before  coupling.The  data  points  are  obtained  from  curve-fitting  the  measured  transmission  spectra  to  the  analytical  model.  a,b,  Comparison 
when  resonators  are  passive  (no  gain  in  pRi:  red  square  and  triangular  symbols)  with  Q-factors  2.9  x  107  and  3.0  x  107,  respectively,  for  pRi  and  pR2,  and 
when  one  resonator  is  active  and  the  other  is  passive  (gain  in  pRp  blue  circular  and  cross  symbols).  c,d,  Effect  of  initial  Q-factor  (loss)  of  pR2  on 
eigenfrequencies.  Two  resonances  with  Q-factors  2.0  x  107  (blue)  and  3.0  x  107(red)  are  chosen  for  pR2.  Shaded  regions  correspond  to  the 
broken-PT-symmetric  region  when  gain  and  loss  are  balanced. The  dashed  and  solid  lines  were  obtained  from  the  theoretical  model  using  experimentally 
obtained  values  of  coupling  strength,  gain  and  the  loss  parameters  (Supplementary  Sections  B2-B4  and  Figs  5  and  6).  The  value  of  Kpj  in  c,d  is  obtained  for 
the  case  when  highQ=  3.0  x  107  resonance  mode  is  chosen  in  pR2. 


broken-symmetry  phase,  due  to  the  stronger  field  localization  in 
the  resonator  with  gain,  as  compared  with  the  unbroken- symmetry 
phase  (Supplementary  Section  7  and  Fig.  9). 

As  a  result  of  stronger  nonlinearity  in  the  broken- symmetry  case, 
the  PT  transition  is  associated  with  a  transition  from  reciprocal 
to  non-reciprocal  behaviour.  When  the  pump  at  port  1  was  OFF 
(liRj  and  |_|R2  are  passive)  and  a  weak  probe  light  was  input 
at  port  1  or  4,  we  observed  resonance  peaks  in  the  forward  or 
backward  transmissions  (Fig.  4a(i),b(i))  with  no  resolvable  mode 
splitting.  When  the  pump  was  set  ON  and  the  gain  and  loss 
were  balanced  so  as  to  operate  in  the  unbroken-PT-symmetric 
region,  transmission  spectra  showed  amplified  signals  with  clearly 
resolved  split  peaks  (Fig.  4a(ii),b(ii)).  Flowever,  when  the  coupling 
strength  was  decreased  so  that  the  system  transited  into  the  broken- 
symmetry  region,  forward  transmission  reduced  to  zero  0 

(Fig.  4a(iii))  but  the  backward  transmission  remained  high 
(Fig.  4b(iii)).  The  transmission  spectra  showed  a  single  resonance 
peak,  as  expected  from  the  theory.  Thus,  in  the  broken- symmetry 
region  the  input  at  port  4  was  transmitted  to  port  1  at  resonance; 
however,  the  input  at  port  1  could  not  be  transmitted  to  port  4,  in 
stark  contrast  with  what  was  observed  for  the  unbroken- symmetry 
region.  This  indicates  non-reciprocal  light  transport  between  ports 
1  and  4. 

Unlike  previous  experiments  demonstrating  non-reciprocal 
transport  in  non-PT  structures34^38  and  asymmetric  behaviour  in 
PT  electronics13,  we  observed  a  complete  absence  of  resonance  peak 
in  the  forward  transmission.  The  advantages  off  he  present  design, 
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which  brings  together  PT- symmetric  concepts  with  nonlinearity- 
induced  non-reciprocal  light  transmission,  over  the  non-PT 
schemes  utilizing  nonlinearity  are:  a  significant  reduction  in  the 
input  power  to  observe  non-reciprocity  (B  1  |JW  in  this  work  versus 
3  W  in  ref.  34,  0.310  W  in  ref.  35  and  B  85  |iW  in  ref.  37);  higher 
contrast;  small  footprint;  and  a  complete  absence  of  the  signal  in 
one  direction  but  resonantly  enhanced  transmission  in  the  other 
direction.  This  is  in  stark  contrast  with  other  non-reciprocal  devices 
where  transmission  suffers  from  high  losses  (see  Supplementary 
Section  Bll  for  a  detailed  comparison  off  his  work  with  other 
works  in  the  literature).  Also,  the  transmitted  signal  here  was  not 
from  spontaneous  emission  off  he  gain  medium.  Without  the  weak 
injected  signal  at  input  port  4,  the  output  at  port  1  was  at  the  noise 
level,  and  no  resonance  peak  was  observed  (inset  of  Fig.  4b  ( ii)  ,b  ( iii) ) . 
Resonance  enhancement  (the  peak)  was  observed  only  when  the 
weak  signal  was  present  (Supplementary  Section  B8  and  Fig.  10).  We 
also  observed  similar  non-reciprocity  between  ports  2  and  3.  These 
results  imply  that  PT-symmetric  WGMRs  can  have  strong  non¬ 
reciprocal  effects  (all-optical  diode  action)  in  the  nonlinear  regime 
with  a  very  low  power  threshold  due  to  significant  enhancement  of 
nonlinearity  in  the  broken- symmetry  phase. 

Our  experiments  extend  PT-symmetric  optics  from  centimetre/ 
metre-scale  structures  to  on-chip  micro-scale  structures  and,  more 
importantly,  from  waveguides  to  microresonators.  This  overcomes 
the  long-standing  open  challenge  of  implementing  a  true  PT- 
symmetric  resonator  system  in  the  optical  regime  and  opens 
the  door  for  exploring  new  functionalities  (for  example, 
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Figure  3  |  Input-output  relation  in  PT-symmetric  WGMRs  and  reciprocity  in  the  linear  regime,  a,  Linear  input-outputrelation  in  the  unbroken-symmetry 
region  (red  square  symbols)  but  nonlinear  in  the  broken-symmetry  region  (blue  circle  symbols).  Transmitted  signal  to  port  1  was  detected  by  a 
photodetectorwhen  the  input  was  at  port  4.  Each  data  pointis  the  average  often  measurements,  and  the  error  bars  represent  the  standard  deviations. 
b,c,  Transmission  spectra  in  the  linear  regime(input  power  of  B  80  nW)  show  reciprocal  lighttransmission  at  forward  (blue  inset)  and  backward  (red 
spectra)  directions  in  both  the  broken-  (b)  and  unbroken-  (c)  symmetry  regions.  The  red  (blue)  spectra  were  normalized  with  the  signal  detected  at  port  3 
(2)  when  the  input  was  at  port  4  (1)  and  there  was  no  coupling  between  the  fibre  tapers  and  the  resonators.  The  slight  difference  in  the  heights  of  the 
resonance  peaks  is  attributed  to  the  laser  power  fluctuations  during  frequency  scanning  and  to  the  thermal  fluctuations  of  the  environment. These  can  be 
minimized,  if  not  completely  eliminated,  by  active  stabilization  and  better  control  of  the  experimental  conditions  (Supplementary  Section  B10  and  Fig.  12). 
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Figure 4  |  Experimentally  observed  unidirectional  transmission  for  PT-symmetric  WG  M  microresonators  in  the  nonlinear  regime.  a,b,  When  both 
resonators  are  passive  (zero  gain),  the  transmission  is  bi-directional  (reciprocal),  and  light  is  transmitted  in  both  forward  (a(i))  and  backward  directions 
(b(i)).  In  the  unbroken-symmetry  region,  where  the  coupling  exceeds  the  critical  value  and  gain  and  loss  are  balanced,  the  transmission  is  still 
bi-directional  (a(ii),  b(ii)).  Mode  splitting  due  to  coupling  is  now  resolved  because  gain  compensates  loss  leading  to  narrower  linewidths.  In  the 
broken-symmetry  region  (a(iii),  b(iii)),  transmission  becomes  unidirectional  (non-reciprocal).  Input  in  the  backward  direction  reaches  the  output  (b(iii)), 
but  input  in  the  forward  direction  does  not  (a(iii)).  This  resembles  the  action  of  a  diodeand  implies  that  an  all-optical  on-chip  diode  with  PT-symmetric 
WG  M  microcavities  operates  in  the  broken-symmetry  region. Inset  in  (b(ii),  b(iii))  shows  the  signal  at  port  1  when  there  is  no  input  signal  at  port  4. 


coherent-perfect-absorption  lasers,  topologically  protected  optical 
diodes,  enhanced  nonlinearities  and  light-matter  interactions) 
that  can  be  achieved  only  with  resonant  structures  and  resonant 
enhancement.  Although  we  used  fibre-taper-coupled  microtoroids, 
our  techniques  can  be  extended  to  other  WGMRs  and  to  photonic 
crystal  cavities.  Similarly,  gain  could  be  provided  by  quantum  dots 
or  other  rare-earth  ions,  and  also  through  nonlinear  processes,  such 
as  Raman  or  parametric  amplification.  Like  any  non-reciprocal 
device  utilizing  resonant  effects,  our  PT-symmetric  all-optical 
diode  is  bandwidth-limited36'38.  However,  by  thermally  tuning 
resonance  wavelengths  and  by  using  active  resonators  doped  with 
multiple  rare-earth  ions,  operation  over  large  wavelength  bands 


should  be  possible.  Coupled  WGMRs  provide  a  comprehensive 
framework  for  understanding  resonance  effects  in  PT-symmetric 
optical  systems  and  could  thereby  aid  in  developing  on-chip 
synthetic  structures  to  harness  the  flow  of  light.  For  example, 
the  electromagnetically  induced  transparency  in  coupled  passive 
resonators  may  benefit  from  PT-symmetric  resonators  through 
lossless  modulation  oft  he  transparency  for  slowing  and  stopping 
of  light.  Similarly,  PT-symmetric  microresonators  can  be  used 
for  studying  nonlinear  Fano  resonances  that  may  give  rise  to 
ultralow-power  and  high-contrast  switching  and  non-reciprocity 
due  to  their  sharp  asymmetric  line  shapes.  Moreover,  there  has  been 
an  emerging  interest  in  exploring  PT  symmetry  in  various  fields, 
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such  as  microlasers42,  sensing46-49,  plasmonics48,50,  optomechanics51 
and  cavity-quantum  electrodynamics52,  where  passive  WGMRs 
have  been  traditionally  used.  This  may  bring  about  new  results  and 
physical  insights  into  these  fields. 
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Parity-time-symmetric  whispering-gallery  microcavities 


Al.  F abrication  of  Erbium-doped  silica  microtoroid.  The  PT-symmetric  optical  resonator 
system  investigated  here  is  formed  by  two  coupled  WGM  microtoroids.  The  passive  (lossy) 
microtoroid  is  made  from  silica1,  whereas  the  active  (with  gain)  resonator  is  made  from 
erbium-ion-doped  silica  film,  which  is  formed  using  sol-gel  synthesis2'3.  The  flow  chart  for 
the  preparation  of  the  sol-gel  silica  film  is  as  follows  (Fig.  SI): 

1)  Sol-gel  precursor  solution  is  prepared  by  mixing  tetraethoxysilane  (TEOS)  in  isopropanol 
alcohol  (IPA),  water  (H20),  and  hydrochloric  acid  (HC1).  The  weight  ratio  of 
TE0S:IPA:H20:HC1  is  0.6:6.5:0.7:6.1. 

2)  Erbium  ions  are  incorporated  by  adding  erbium  nitrate  Er(N02)3  to  the  precursor  solution 
at  desired  concentration. 

3)  The  solution  is  stirred  for  3  hours  at  70°C  during  which  hydrolysis  and  condensation 
occur. 

4)  The  solution  is  aged  for  24  hours  at  room  temperature. 

5)  The  aged  solution  is  then  spin-coated  on  the  silicon  wafer  to  form  a  uniform  layer. 

6)  The  film  is  heat  treated  at  1000°C  for  3  hours  to  form  an  erbium-doped  silica  glass. 

The  process  is  repeated  multiple  times  until  the  desired  silica-film  thickness  is  obtained.  Note 
that  without  step  2),  the  resultant  film  would  be  a  pure  sol-gel  silica  film. 


A2.  F abrication  of  edge  microtoroids.  To  investigate  the  PT-symmetric  coupled-resonator 
system,  we  must  be  able  to  tune  the  coupling  strength  and/or  the  gain-to-loss  ratio.  The 
optical  gain  in  the  active  resonator,  which  enables  us  to  tune  the  gain-to-loss  ratio,  is 
provided  by  optically  pumping  the  erbium  ions  with  a  pump  laser45.  The  coupling  strength  is 
tuned  by  varying  the  distance  between  the  two  resonators6.  This  can  be  achieved  by 
fabricating  each  of  the  resonators  at  the  edge  of  a  different  chip  and  by  controlling  the 
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FIG.  SI.  Preparation  of  erbium-doped  silica  film  for  fabrication  of  active  WGM 
resonators.  Any  rare-earth  ion,  such  as  ytterbium,  neodymium,  and  thulium,  can  be 
incorporated  into  the  silica  matrix  using  this  process  by  replacing  the  lir( NO  3 f  with  a 
soluble  compound  containing  the  relevant  ion.  Without  the  additional  step  for  rare-earth-ion 
or  other  gain  material,  the  resultant  film  would  be  a  plain  sol-gel  silica  film. 
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FIG.  S2.  Fabrication  process  of  edge -microtoroids  on  silicon  wafer.  Top  (upper  panel )  and 
the  side  (lower  panel )  views  of  the  resultant  structure  during  the  fabrication  are  shown.  The 
fabrication  steps  1-5,  9  and  10  are  the  steps  for  standard  silica  microtoroid  fabrication  on 
silicon  wafer.  The  steps  in  the  dotted  box  are  necessary  for  obtaining  the  edge  microtoroids, 
which  stand  on  silicon  pillars  and  hang  over  the  silicon  wafer.  A  passive  silica  edge- 
microtoroid  is  fabricated  starting  in  step  1  with  a  thermal  silica  or  sol-gel  silica  film  on 
silicon  wafer.  Starting  in  step  1  with  rare-earth-ion  doped  sol-gel  silica  film  on  silicon  wafer, 
an  edge  active-microtoroid  is  fabricated .  Si:  silicon;  .S'/CT:  silica;  PR:  photoresist;  HF: 
hydrofluoric  acid;  Xeh'2:  xenon  difluoride;  CO2:  carbon  dioxide  (see  text  for  details). 


separation  between  chips  by  using  nanopositioning  systems  on  which  the  chips  are  placed. 

To  fabricate  the  microtoroids  at  the  edges  of  the  chips,  we  modify  slightly  the  original  recipe 
for  fabricating  on-chip  microtoroid  resonators1'7.  The  process,  which  is  illustrated  in  Fig.S2, 
is  the  same  for  both  passive  and  active  resonators.  It  begins  with  silica  free  of  gain  dopants 
on  a  silicon  wafer  for  the  passive  resonator  and  with  erbium-doped  silica  film  on  a  silicon 
wafer  for  the  active  one.  We  fabricate  the  edge-toroids  as  follows: 
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1)  A  photoresist  (PR)  layer  is  spin-coated  over  plain  silica  (for  the  passive  resonator)  or 
erbium-doped  sol-gel  silica  (for  the  active  resonator). 

2)  Using  UV-photolithography  circular  disks  are  patterned  on  the  silica  film. 

3)  PR  is  then  developed,  forming  PR  disks. 

4)  With  hydrofluoric  (HF)  acid  as  the  etchant,  silica  that  is  not  covered  with  the  PR  is 
removed  in  order  to  form  PR-coated  silica  disks  on  silicon  wafer. 

5)  PR  is  then  removed  by  washing  the  wafer  with  acetone ,  uncovering  the  silica  disks. 

6)  and  7)  A  new  layer  of  PR  is  spun  coated  on  the  wafer  and  then  the  chip  is  exposed  to  XeF2 
gas,  which  isotropically  etches  silicon.  The  PR  layer  forms  a  protective  layer,  so  XeF2  does 
not  etch  the  structure  from  the  top.  Etching  only  proceeds  in  a  direction  parallel  to  the 
surface . 

8)  The  PR  is  washed  away  with  acetone.  Steps  6)-8)  are  repeated  until  the  desired  over-hang 
disk  structure  is  formed. 

9)  The  wafer  is  immersed  in  XeF2  gas  once  more  to  etch  silicon  from  the  top  and  sides  in 
order  to  form  the  pillar  structure,  i.e.,  silica  disks  over  silicon  pillars. 

10)  Finally,  CCF  reflow  heats  and  melts  the  silica  disks,  transforming  them  into  silica 
microtoroids.  The  resulting  structures  are  microtoroids  at  the  edge  of  a  silicon  wafer  with 
their  pillars  on  the  silicon  substrate  but  with  a  portion  of  the  silica  torus  extending  beyond  the 
wafer. 


A3.  Fabrication  of  fiber  tapers  for  add-drop  filter  (ADF)  configuration.  A  pair  of  tapered 
optical  fibers  is  used  to  couple  light  into  and  out  of  the  WGM  resonators.  The  coupled 
resonator  system  is  placed  between  two  fiber  tapers,  similar  to  the  standard  add-drop-filter 
(ADF)  configuration  (Fig.  S3  and  Fig.  S4B).  The  fiber  tapers  should  be  parallel  to  each  other 
and  have  similar  waist  size  and  mode  profiles  to  achieve  the  best  coupling.  Their  separation 
must  be  large  enough  to  place  the  coupled  resonators  between  them.  To  achieve  these,  the 
fiber  tapers  were  fabricated  by  heating  and  stretching  a  pair  of  single  mode  silica  fibers 
simultaneously  over  a  hydrogen  flame.  The  result  was  two  fiber  tapers  with  a  very  small 
height  difference  and  ~90%  transmission.  We  used  a  fiber-tip  controlled  with  a 
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nanopositioner  to  push  one  of  the  fiber  tapers  away  from  the  other  until  their  separation  was 
roughly  equal  to  the  sum  of  the  diameters  of  the  two  resonators.  The  separation  was  adjusted 
by  pushing  the  tapers  towards  or  away  from  each  other  as  well  as  to  and  from  the  resonators 
until  the  desired  coupling  conditions  were  achieved8.  An  additional  nanopositioning  stage 
was  used  to  adjust  the  distance  between  the  resonators. 


B.  Experimental  setup  and  characterization  of  the  PT-symmetric  microcavities 


Bl.  Experimental  setup.  Our  experimental  setup  consists  of  seven  ingredients  (Fig.  S3). 

1)  The  coupled  system  of  a  passive  and  an  active  microtoroid.  The  resonators  are 
directly  coupled  to  each  other  in  the  1550  nm  band  via  their  evanescent  fields.  There  is  no 
direct  coupling  between  the  resonators  in  the  1460  nm  band. 

2)  A  pair  of  tapered  optical  fibers  to  couple  light  into  and  out  of  the  WGM  modes  of 
the  resonator.  Each  tapered  fiber  is  coupled  to  one  resonator. 


1450nm 

Tunable 

Laser 


Computer 


FIG.  S3.  Experimental  setup  used  for  the  study  of  PT-symmetric  whispering  gallery  mode 
(WGM)  microcavities.  uRl :  active  microtoroid;  /uR2:  passive  microtoroid;  PD: 
photodetector;  WDM:  wavelength  division  multiplexer;  FC:  fiber  connector ;  TEC:  thermo¬ 
electric  cooler;  OSA:  optical  spectrum  analyzer. 
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FIG.  S4.  Illustrations  of  the  setups  used  in  the  experiments.  (A)  Setup  used  to  show  the  PT 
phase  transition  presented  in  Fig.  2  in  the  main  text.  [A(i)  and  A(ii)]  correspond  to  two  cases 
where  the  pump  is  input  in  the  same  (or  opposite)  direction  as  the  probe  signal,  respectively . 
The  results  remain  the  same  regardless  of  the  pump  direction.  (B)  Setup  used  for  the 
demonstration  of  nonreciprocal  light  transmission  as  shown  in  Figs.  3  and  4  of  the  main  text. 
[B(i)  and  B(ii)J  respectively  correspond  to  forward  and  backward  transmission.  pRp  active 
microtoroid;  11R2:  passive  microtoroid;  PD:  photodetector .  Signal  probe  and  pump  fields  are 
denoted  with  different  colors. 


3)  Two  external-cavity  tunable  laser  sources.  The  first  is  used  as  the  probe  (signal)  in 
the  1550  nm  band;  the  second  is  in  the  1460  nm  band  and  is  used  to  pump  the  erbium  ions  in 
the  active  microtoroid. 

4)  Nanopositioning  systems  to  tune  and  adjust  the  coupling  strengths  between  the 
resonators  and  fiber  tapers  as  well  as  between  the  directly-coupled  active  and  passive 
resonators. 

5)  Thermoelectric  cooler  (TEC)  used  to  tune  the  resonance  of  the  passive  resonator 
via  the  thermo-optic  effect.  Initially,  the  resonances  of  the  two  resonators  in  the  1550  nm 
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band  are  spectrally  separated  from  each  other.  Using  the  TEC,  we  tune  one  of  the  resonance 
lines  of  the  passive  resonator  to  overlap  with  a  resonance  line  of  the  active  resonator  in  the 
same  band  [6],  The  resonators  do  not  have  overlapping  resonance  lines  in  the  1460  nm  band, 
so  the  light  from  the  pump  laser  is  coupled  only  to  the  active  resonator. 

6)  Photodetectors  (PD)  connected  to  an  oscilloscope  to  monitor  the  transmission 
spectra  at  the  outputs.  Pump  and  probe  light  are  separated  using  wavelength  division 
multiplexers.  An  optical  spectrum  analyzer  is  used  to  characterize  the  optical  gain  and  lasing 
from  erbium  ions. 

7)  A  computer  to  process  the  transmission  spectra  and  extract  information  on  mode 
splitting  (difference  between  the  real  parts  of  the  eigenfrequencies  of  the  coupled  system)  and 
the  linewidth  difference  (difference  between  the  imaginary  parts  of  the  eigenfrequencies  of 
the  coupled  system)  of  the  resonance  lines. 

We  performed  two  sets  of  experiments  using  the  setup  in  Fig.  S3  in  different 
subconfigurations  by  proper  connections  of  the  input  and  output  ports: 

a)  By  moving  one  of  the  fiber  tapers  far  from  the  coupled  system,  we  determined  the 
unbroken-  and  broken-symmetry  regions  as  a  function  of  the  coupling  strength  between  the 
resonators  and  as  a  function  of  the  gain-to-loss  ratio  (determined  by  the  pump  power).  The 
results  are  shown  in  Fig.  2  of  the  main  text.  Figure  S4A  illustrates  the  scheme  used  here. 
Experiments  were  performed  when  the  pump  and  probe  (signal)  are  input  in  the  same 
direction  [Fig.  S4A(i)]  and  when  they  are  in  opposite  directions  [Fig.  S4A(ii)]. 

b)  By  using  the  ADF  configuration,  which  is  essentially  a  2-by-2  optical  device,  we 
studied  unidirectional  light  transmission  (nonreciprocal  optical  transport).  Figure  S4B  depicts 
the  illustration  of  the  scheme  used  here.  The  pump  laser  was  input  at  port  1  and  the  probe 
(signal)  was  input  at  port  1  or  at  port  4.  Transmission  of  the  signal  from  port  1  to  port  4  is  the 
forward  transmission  whereas  transmission  from  port  4  to  port  1  is  the  backward 
transmission. 


B2.  Determining  the  coupling  strength  k  between  the  resonators.  The  coupling  strength 
between  the  resonators  is  a  function  of  the  distance  between  them.  In  our  experiments,  we 
tune  the  distance  between  the  resonators  finely  using  a  nanopositioning  system  with  a  step 
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resolution  of  100  nm.  In  the  following,  we  will  introduce  the  theoretical  framework  for  the 
determination  and  calibration  of  the  coupling  strength  from  the  distance  between  the 
resonators.  We  performed  the  experiment  using  the  configuration  given  in  Fig.  S4A(i)  in 
which  one  of  the  resonators  is  directly  coupled  to  an  optical  fiber  waveguide.  Defining  the 
intracavity  mode  fields  of  the  resonators  as  ak=12  for  the  first  and  second  resonators  with 

resonance  frequencies  cok=l  2  ,  the  coupling  strength  between  the  resonators  as  k  ,  and  the 
input  field  as  ain ,  we  can  write  the  following  rate  equations  for  the  coupled-resonators  system 


-77  =  -to i«,  -  «i  -  iKa2  -  ^Ka,n 

dt  2  (S.l) 

da ,  .  y, 

— -  =  -ico^a,  — -a,  -  iKa . 
dt  2  -  2  -  1 


together  with  the  input-output  relations  aout  =  am  +  .  Here  yk  denotes  the  loss  or  gain  of 

the  resonators  (loss:  77  >0;  gain:  yk  <0),  and  yc  >0  corresponds  to  the  coupling  loss  between 
the  first  resonator  and  the  fiber  taper  waveguide.  The  rate  equations  take  the  form 


=  (A  -  -  ix A  -  n , 

^  2  (S.2) 

^  =  (iA2  -  A)  a,  -  i/cAt 
dt  2 

after  substituting  ak  =  \e~"“  =  1,2)  and  +  ^Le~i01>  ■  Similarly,  for  the 

dt  dt 

input-output  relation  we  have  Aout  -  Ain  +  ^/TA  •  Here,  A;h  2  =co-cok  is  the  detuning  between 
the  resonance  frequencies  and  the  frequency  of  the  input  laser  light.  The  coupling  of  these 
two  resonators  creates  two  supermodes  A+  =(4,  +  AAl-Jl  and  h  =  ( h,  -  A, )  /  V2  with  the 
eigenfrequencies  co+  and  co  given  as 


1  . 

co,  -  —  1 
±  2 


-i  {ay  +a>2 ) 


Yx+Yc+Yi 


-co2)~ 


YAYc-Yi 


(S.3) 


Since  in  the  experiments  we  tune  the  resonances  of  the  resonators  to  be  degenerate 
oy  =  co2  =  cot] ,  the  eigenfrequencies  can  be  re-written  as 


co. 


co, 


h-fa+Yc+Yi) 


±^6A-(7,  +YC-Y2)1 


(S.4) 
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where  the  expression  in  the  square -root  quantifies  the  effect  of  the  coupling  and  the  interplay 
between  the  coupling  strength  and  the  loss/gain  in  the  resonators.  Here  we  define  the 
difference  between  the  eigenfrequencies  as  the  spectral  distance 


(S.5) 


There  are  four  different  regimes  which  are  characterized  by  the  interplay  between  the 
coupling  strength  and  the  total  loss  of  the  system,  (i)  No  coupling  between  the  resonators 
k  =  0.  The  resonance  frequencies  of  individual  resonators  are  equal  and  have  the  value  a>0 , 

but  their  losses  are  different  and  quantified  as  yx  +  y  and  y2  for  the  first  and  second 

resonators,  respectively,  (ii)  Strong  coupling  regime  where  16/c2  >  (yx  +  yc  - y2f  is  satisfied. 

In  this  case,  the  coupling  between  the  resonators  creates  two  supermodes,  whose  frequencies 
are  up-  and  down-shifted  by  SI 2  with  respect  to  the  initial  resonance  frequency  a>0  of 

uncoupled  individual  resonators,  implying  that  supermodes  are  identified  by  the  difference  of 
their  resonance  frequencies.  The  transmission  spectrum  should  exhibit  two  resonances  (mode 
splitting)  whose  frequency  difference  (amount  of  splitting)  increases  with  increasing  coupling 
strength  k  .  Coupling  does  not  induce  extra  loss  and  the  linewidths  of  the  supermodes  are 
quantified  by  the  sum  of  the  losses  of  individual  resonators,  (iii)  Critical  coupling  regime 
where  16/c2  =  (yt  +  yc  -  y2f  .  In  this  regime  S'  =  0  and  the  resonance  frequencies  of  the 

supermodes  coalesce,  (iv)  Weak  coupling  regime  where  16/c2  < (yx  +  yc  - y2)2  is  satisfied. 

Here,  coupling  induces  extra  loss  to  the  system,  and  the  two  supermodes  are  identified  by  the 
difference  in  their  linewidths. 

We  measure  yk  >  0  and  yc  >  0  in  the  experiments  from  the  resonances  in  the  transmission 
spectra  (see  Fig.  Id  of  the  main  text)  of  each  resonator  separately  (i.e,  without  coupling 
between  them).  Note  that  yk  is  related  with  all  the  losses  (scattering  loss,  bending  loss  and 
material  loss)  in  the  resonators  except  for  the  coupling  losses.  Thus  we  need  to  minimize  the 
effect  of  the  coupling  losses  on  our  transmission  spectra  measurements  to  derive  yk .  This  is 

possible  when  the  fiber-resonator  coupling  is  set  to  deep-undercoupling  regime  where 
coupling  loss  is  much  smaller  than  the  other  losses  in  the  system.  In  the  deep -undercoupling 
regime,  we  obtain  the  transmission  spectra  for  each  of  the  resonators  separately  and  then  find 
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the  linewidths  A cok  and  the  resonance  frequency  cok  of  the  mode  of  interest  by  Lorentzian 
curve  fitting  to  the  resonances  obtained  in  the  transmission  spectra.  Then  we  calculate  the 
quality  factor  Qk=I  ,  =  o\  /  A  ay  =  ay  /  yk  from  which  the  loss  parameter  yk  is  calculated.  In 

the  configuration  shown  in  Fig.  S4A(i)  and  in  the  model  expressed  in  Eq.  S.l,  the  first 
resonator  is  coupled  to  a  fiber  taper;  therefore,  the  resonator-fiber  coupling  loss  should  also 
be  estimated  in  the  experimental  settings.  In  order  to  do  this,  we  move  the  fiber  taper  from 
the  deep-undercoupling  regime  closer  to  the  critical  coupling  condition.  This  naturally 
induces  loss  and  leads  to  the  broadening  of  the  linewidth  of  the  resonance  mode  observed  in 
the  transmission  spectrum  of  the  first  resonator.  Under  this  condition,  the  measured  quality 
factor  is  generally  referred  to  as  the  loaded  quality  factor  and  expressed  as  Ql'aded  =  Qt  +  Q~c  1 
where  Q  is  the  fiber-resonator  coupling  quality  factor  and  denotes  the  coupling  loss.  In  other 
words,  loss  parameters  of  the  system  satisfy  yloaded  =  yx  +  yc  ■  Plugging  in  this  expression  the 
value  of  Q,  (or  yl )  estimated  in  the  deep-undercoupling  regime  for  the  first  resonator  and  the 
newly  estimated  Qloaded  (or  yloaded ),  we  find  Q  (or  yc )  used  in  our  experiments.  Thus  in  the 


FIG.  S5.  Relation  between  the  coupling  strength  k  and  the  distance  between  the 
microtoroid  resonators.  The  zero  distance  here  corresponds  to  the  case  when  the  resonators 
are  close  to  each  other  within  the  lOOnm  step-resolution  of  the  nanopositioning  system. 
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expressions  in  Eqs.  S.1-S.4,  the  only  unknown  parameter  is  the  coupling  strength  k  between 
the  resonators.  Note  that  in  the  case  of  optical  pumping  of  the  active  resonator,  photons 
emitted  from  the  gain  medium  provide  gain  to  the  resonance  mode  of  the  first  resonator, 
leading  an  increase  in  Qx  and  decrease  in  yx  as  seen  in  Fig.le  of  the  main  text. 

In  the  experiments,  we  monitored  the  change  in  the  mode  splitting  (difference  between  the 
resonance  frequencies  of  the  supermodes)  as  the  distance  between  the  resonators  is  varied. 
From  the  experimentally  obtained  splitting,  we  estimated  the  value  of  k  using  Eq.  S.4.  The 
resultant  relation  between  k  and  the  distance  between  the  resonators  is  given  in  Fig.S5  where 
we  see  that  the  coupling  strength  exponentially  decreases  with  increasing  distance  between 
the  resonators.  This  result  agrees  well  with  previous  reports  in  the  literature6'910.  A  curve 
fitting  method  to  estimate  the  coupling  strength  k  will  be  explained  in  the  next  section. 


B3.  Estimating  eigenfrequencies  of  the  supermodes  from  the  experimentally  obtained 
transmission  spectra.  Using  the  configuration  given  in  Fig.  S4A(i),  we  probed  the 
transmission  spectra  of  PT-symmetric  coupled  resonators  system  at  different  coupling 
conditions,  and  estimated  the  eigenfrequencies  of  the  supermodes  from  the  measured 
transmission  spectra.  This  is  done  by  curve  fitting  an  analytical  expression  obtained  for 
transmission  using  coupled-mode  theory  to  the  experimentally  obtained  transmission  spectra. 
Starting  with  the  expressions  given  in  Eqs.  S.1-S.2  and  the  input-output  relation  for  the 
system,  we  find  the  steady-state  normalized  transmission  T  =  \Aout  I Ainf  as 


_ 2A(r2-?'2A) _ 

4k2  +  (y,  +yc -i2A)(y2  -i'2A) 


(S.5) 


where  yx ,  y2  and  yc  are  all  known  from  the  experiments  as  explained  above ,  and  we  used 
a>x=Q)2=c y0  and  A  =  co-ay  because  in  the  experiments  the  resonance  frequencies  of  the 

resonators  are  set  to  be  equal  using  thermal  tuning  via  the  thermo-optic  effect.  We  used  the 
expression  in  Eq.  S.5  as  the  curve  fitting  function  to  the  experimentally-obtained 
transmission  spectra  under  different  coupling  strengths  k  (i.e.,  distance  between  the 
resonators).  Note  that  for  the  ideal  PT-symmetry  case  in  which  gain  in  the  active  resonator 
equals  to  loss  in  the  passive  one  (i.e.,  y1+yc<  0  ,  y2>  0  and  yt  +  yc:  =  -y2  ),  Eq.  S.5  becomes 
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2rc{r2-'2A) 

4  k1  -  {y\  +  4A2) 


(S.6) 


When  there  is  no  gain  in  the  system  (coupled  passive  resonators),  all  the  parameters  except 
for  are  known.  Therefore,  curve  fitting  the  expression  in  Eq.  S.5  to  the  experimentally - 
obtained  transmission  spectrum  with  k  as  the  free  parameter  provides  a  value  for  k  which  is 
then  used  in  Eq.  S.4  to  obtain  the  eigenfrequencies  of  the  two  supermodes.  When  the  gain  is 
present  in  the  system,  we  used  k  and  y{  as  the  free  parameters  for  curve  fitting.  The 

estimated  values  of  k and  ;/,  are  used  in  Eq.  S.4  to  obtain  the  supermode  eigenfrequencies. 


In  Fig.  S6,  we  provide  the  experimentally  obtained  transmission  spectra  together  with  the 
fitted  curves  using  Eq.  S.5  and  the  eigenfrequencies  estimated  via  Eq.  S.4.  Clearly,  even  if 
the  resonances  coalesce  as  shown  in  Fig.  S6A,  the  theoretical  fit  provides  two  different 
eigenfrequencies  distinguished  by  the  difference  in  their  line  widths.  The  confidence  of  the 
curve  fitting  was  always  larger  than  0.95. 
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FIG.  S6.  Theoretical  curve  fitting  to  extract  the  eigenfrequencies  in  the  broken-PT- 
symmetric  phase.  We  fitted  the  transmission  function  obtained  from  the  theoretical  model 
and  given  in  Eq.  S.5  to  the  transmission  spectra  obtained  in  the  experiments .  Estimated 
coupling  strengths  tc  and  gain-loss  parameter  yx  are  (A)  k  =  295  MHz  and 
;/,  =-13.84  MHz.  ,(B)  ft-  =  4.16  MHz  and  y  =-13.07  MHz. ,  and(C)  k  =  9.75  MHz  and 
yl  =-14.55  MHz  •  The  detunings  of  the  calculated  eigenfrequencies  from  the  initial  resonance 
frequencies  are  given  in  the  inset  and  the  units  are  in  MHz.  Blue  denotes  the  experimental 
data  and  red  denotes  the  best  fit. 
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B4.  Threshold  coupling  strength  kpt  .  In  Fig.  2  of  the  main  text  and  Fig.  S7  of  this 
Supplement,  we  plotted  the  real  and  imaginary  parts  of  the  eigenfrequencies  as  a  function  of 
the  parameter  k!  kpt  where  at  is  the  coupling  strength  as  discussed  above.  Here  kpt  is  the 

threshold  coupling  strength  at  which  the  real  parts  of  the  eigenfrequencies  of  the  supermodes 
coalesce  for  the  first  time  and  become  equal  for  a  fixed  gain  and  loss  ratio  as  the  coupling 
strength  a:  (i.e.,  distance  between  the  resonators)  is  decreased.  Then  the  value  of  k  satisfying 
16 k1  =  (;/,  +  yc  -y2  )2  is  the  threshold  coupling  strength  kpt  ,  which  is  found  as 


1 , 

kpt  =-\ri+rc-r2 


(S.7) 


It  is  clear  that  kpt  depends  on  the  gain-loss  ratio  of  the  coupled  system.  For  the  ideal  case 
when  gain  and  loss  are  balanced  (i.e.,  yl+yc=-y2),  the  threshold  coupling  constant  kpt 
becomes 


1 

KPT 


(S.8) 


If  one  can  measure  y1 ,  yc  and  y2  precisely  then  the  kpt  can  be  calculated  using  Eq.  S.7. 

However,  in  an  experiment,  especially  when  there  is  gain  in  the  system,  one  does  not  know 
the  exact  value  of  yx  and  may  not  achieve  exactly  balanced  loss  and  gain.  Moreover,  in 

scheme  like  ours,  the  coupling  strength  that  can  be  probed  is  limited  by  the  step  resolution  of 
100  nm  of  the  nanopositioning  system.  Therefore,  satisfying  the  above  expression  exactly  is 
difficult. 

In  the  experiments,  we  first  bring  the  loss  and  the  gain  as  close  to  each  other  as  possible,  and 
then  monitor  the  transmission  spectra  as  we  decrease  the  coupling  strength  k  by  increasing 
the  distance  between  the  resonators.  At  each  value  of  distance  (i.e.,  our  nanopositioning 
system  has  a  step  resolution  of  100  nm),  we  obtain  the  transmission  spectrum  and  then 
estimate  the  coupling  strength  k  and  yl  via  the  curve  fitting  as  described  in  section  B3.  Then 

using  the  estimated  value  of  k  and  yx  ,  we  calculate  the  real  and  imaginary  parts  of  the 

eigenfrequencies.  The  value  of  k  at  which  the  real  parts  of  the  eigenfrequencies  become  the 
same  for  the  first  time  as  k  is  decreased  (i.e.,  distance  between  the  resonators  is  increased)  is 
the  value  of  kpt  .  Due  to  the  limited  resolution  of  our  nanopositioning  system,  in  the 

experiments  we  may  miss  the  exact  point  where  16/c2  =(yl  +/c  —  is  satisfied. 


13 


NATURE  PHYSICS  |  www.nature.com/naturephysics 


©  2014  Macmillan  Publishers  Limited.  All  rights  reserved. 


SUPPLEMENTARY  INFORMATION 


DOI:  10.1038/NPHYS2927 


Therefore,  we  need  to  estimate  kpt  from  the  available  information.  In  the  experiments,  we 
found  out  that  the  value  of  y]  is  almost  constant  as  it  should  be  because  we  keep  the  gain 
condition  constant  by  fixing  the  pump  power  of  the  gain  medium.  Slight  fluctuations  yx 

come  from  fluctuations  of  the  pump  power  and  frequency.  From  these  experiments,  we 
assigned  an  average  value  toyx .  Then  we  used  this  value  in  Eq.S.7  to  calculate  kpt  .  We 

checked  the  validity  of  the  calculated  kpt  by  confirming  that  it  falls  between  two 

experimentally  accessible  k  values  kx  and  k2  which  satisfy  16xf  >  (y[  +yc  -y2)2  and 

16/c22  <  (yx  +yc -y2)2 ,  (i.e.,  kx  >  kpt  >  k2). 

In  Figs.  2c  and  2d,  we  compare  the  evolution  of  the  real  and  imaginary  parts  of  the 
eigenfrequencies  of  the  supermodes  for  two  different  values  of  the  loss  in  the  second 
resonator.  The  value  of  kpt  used  in  these  figures  was  obtained  for  the  case  of  high  quality 

factor  Q=3.0xl07  resonance  mode  (resonance  mode  with  lower  loss)  of  the  second  resonator 
]iR2  in  order  to  visually  compare  the  symmetry  breaking  points  of  these  two  different  cases. 
Therefore,  the  symmetry  breaking  point  for  the  case  of  lower  loss  appears  at  k/kpt=  1 

whereas  that  for  the  case  with  higher  loss  appears  at  k!  kpt  >  1 .  This  confirms  that  the  higher 

the  loss  in  the  passive  resonator  the  higher  the  threshold  for  symmetry  breaking  point,  as 
quantified  in  Eqs.  S.7  and  S.8. 


B5.  Typical  transmission  spectra  in  the  broken-  and  unbroken-symmetry  regions.  Using 
the  configuration  given  in  Fig.  S4A(i),  we  probed  the  transmission  spectra  of  the  coupled- 
resonators  system  as  a  function  of  the  coupling  strength  after  setting  the  gain  and  loss  as 
close  to  each  other  as  possible.  The  results  showing  the  PT-symmetry  breaking  are  provided 
in  Fig.  2  of  the  main  text.  In  Fig,  S7,  we  provide  typical  transmission  spectra  obtained  in  our 
experiments  for  the  broken-  and  unbroken-symmetry  phases  shown  in  Fig.2  of  the  main  text. 


B6.  Co-  and  counter-propagating  pump  and  probe  light.  In  order  to  check  the  effect  of  the 
direction  of  the  pump  on  the  observed  PT-symmetric  phenomena,  we  performed  experiments 
using  the  configuration  in  Fig.S4A(i)  and  (ii).  Typical  experimentally-obtained  transmission 
spectra  in  the  broken-  and  unbroken-PT  symmetry  phases  are  given  in  Fig.  S8.  It  is  clear  that 
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regardless  of  the  propagation  direction  of  the  pump  laser  with  respect  to  the  probe,  the  same 
behavior  is  observed.  This  can  be  explained  as  follows.  WGM  resonators  support  counter- 
propagating  frequency-degenerate  modes,  so  light  emitted  from  excited  erbium  ions  during 
their  transition  to  the  ground  state  can  couple  to  modes  propagating  in  clockwise  and  counter¬ 
clockwise  directions  with  equal  probability.  Thus,  regardless  of  the  propagation  direction  of 
the  pump,  the  weak  probe  signal  is  always  amplified. 


Fig.S7.  Experimentally  obtained  transmission  spectra  in  broken-PT-  and  unbroke n-PT- 
symmetric  regions.  The  configuration  shown  in  Fig.  S4A(i)  is  used.  (A)  Mode-splitting, 
corresponding  to  the  difference  between  the  real  parts  of  the  eigenfrequencies  of  the 
supermodes  of  the  coupled  WGM-resonator  system,  as  a  function  of  the  coupling  strength. 
The  coupling  conditions  labeled  with  square  marks  and  letters  (B)-(E)  are  the  conditions  for 
which  transmission  spectra  are  given  in  panels  (B)-(E).  The  straight  red  lines  in  (A)  depict 
the  case  in  which  both  resonators  are  passive;  the  dotted  blue  lines  depict  the  case  in  which 
one  resonator  is  active  and  the  other  is  passive  with  balanced  loss  and  gain.  (B)  Both 
resonators  are  passive  (no  gain)  and  the  coupling  is  small.  Splitting  is  barely  seen.  (C)Both 
resonators  are  passive  and  the  coupling  is  strong.  Mode -splitting  is  clearly  seen.  (D) 
Coupled  passive  and  active  resonators  in  the  unbroken- symmetry  region  .  Split  resonance 
peaks  are  clearly  seen.  (E)  Coupled  passive  and  active  resonators  in  the  broken- symmetry 
region,  showing  the  coalescence  of  the  supermodes. 
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Fig.S8.  Experimentally  obtained  transmission  spectra  in  broken-  and  unbroken-PT- 
symmetric  regions  with  different  pump  and  probe  directions.  In  these  experiments  the 
configurations  in  Fig.  S4A  are  used.  (A)  and  (B)are  obtained  in  the  unbroken- symmetry 
region  using  Fig.  S4A(i)  and  Fig.  S4A(ii),  respectively.  (C)  and  (D)are  obtained  in  the 
broken-symmetry  region  using  Fig.  S4A(i)  and  Fig.  S4A(ii),  respectively.  The  transmission 
spectra  are  the  same  regardless  of  the  directions  of  the  pump  and  probe  light  fields. 


B7.  Localization  of  field  in  the  active  resonator  in  the  broken  PT -symmetry  phase.  Strong 
nonreciprocity  observed  in  our  PT-symmetric  optical  resonators  is  due  to  the  significant 
enhancement  of  nonlinearity  in  the  broken  symmetry  phase.  This  enhancement,  on  the  other 
hand,  originated  from  the  strong  field-localization  in  the  broken  symmetry  phase.  Field 
localization  here  means  that  regardless  of  which  port  is  used  as  the  input  port,  the  field  is 
always  localized  in  the  active  resonator  (i.e.,  resonator  with  gain  and  with  less  loss). 
Therefore,  the  signal  at  the  output  port  of  the  fiber  taper  coupled  to  this  active  resonator 
shows  a  strong  resonance  peak  whereas  the  signal  at  the  output  port  of  the  fiber  taper  coupled 
to  the  passive  resonator  (i.e.,  resonator  without  gain)  is  minimized,  if  not  completely 
eliminated.  This  is  true  regardless  of  whether  the  input  is  at  the  fiber  taper  waveguide 
coupled  directly  to  the  passive  or  the  active  resonator.  The  results  of  our  experiments  are 
depicted  in  Fig.  S9  where  we  see  that  only  when  the  PT-symmetry  is  broken,  the  field  is 
localized  in  the  active  resonator  and  thus  the  signal  at  the  output  port  of  the  fiber  coupled  to  it 
shows  a  strong  resonance  peak  whereas  there  is  a  complete  absence  of  resonance  peaks  at  the 
output  of  the  fiber  coupled  to  the  passive  resonator.  In  the  unbroken  phase,  both  outputs  show 
resonance  peaks  regardless  of  the  input  port. 

When  there  is  no  gain  in  the  system  (both  resonators  are  passive),  the  output  of  the  fiber 
through  which  the  field  is  input  shows  a  resonance  dip  and  the  output  of  the  other  fiber  shows 
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a  resonance  peak  [Fig.  S9A(a)  and  S9B(a)  |.  Resonance  dip  at  port  2  (4)  for  the  input  at  port  1 
(3)  is  due  to  destructive  interference  between  the  light  transmitted  directly  to  the  output  and 
the  portion  coupled  to  the  resonator  and  then  coupled  back  to  the  fiber.  A  portion  of  the  light 
coupled  to  the  resonator  from  the  input  fiber  then  couples  to  the  other  resonator  and  leaks  to 
the  output  port  of  the  other  fiber,  leading  to  resonance  peaks. 
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Fig.S9.  Localization  of  the  optical  field  in  the  active  resonator  in  the  broken-PT  symmetry 
phase.  Fiber  taper  waveguide  with  ports  1  and  2  is  coupled  to  the  active  resonator  whereas 
the  one  with  ports  3  and  4  is  coupled  to  the  passive  resonators.  (A)  Input  is  at  port  1  (active 
resonator  side),  and  outputs  are  measured  at  ports  2  and  4.  (B)  Input  is  at  port  3  (passive 
resonator  side),  and  outputs  are  measured  at  ports  2  and  4.  [A(a)]  and  [B(a)]  denote  the  case 
when  there  is  no  gain  in  the  system,  [A(b)]  and  [B(b)]  correspond  to  the  case  where  the  gain 
and  loss  of  the  system  is  balanced  and  the  system  is  in  unbroken-PT-symmetric  phase 
because  coupling  strength  is  large,  and  finally  [A(c)]  and  [B(c)]  correspond  to  the  case 
where  the  gain  and  loss  of  the  system  is  balanced  and  the  system  is  in  the  broken  -PT- 
symmetric  phase  due  to  the  weaker  coupling  strength  between  the  resonators.  Blue  curves 
denote  the  output  at  port  2  and  red  curves  denote  the  output  at  port  4.  It  is  clear  that  in  the 
broken  PT-symmetric  phase,  the  field  is  always  localized  in  the  resonator  with  gain  or  with 
low  loss,  regardless  of  whether  the  input  is  at  the  passive  resonator  side  or  at  the  active 
resonator  side. 
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Fig.SlO.  Experimentally  obtained  spectra  showing  the  effect  of  the  weak  signal  at  the  input 
port  on  the  output  transmission  spectra.  This  is  an  enlarged  version  of  the  spectra  shown  in 
Figs.  4[b(ii)  and  b(iii)]  for  the  backward  transmission  as  depicted  in  Fig.  S4B(ii).  For  both 
the  unbroken-  and  broken-  symmetry  cases,  the  output  port  1  remains  at  the  noise  level 
without  any  observable  resonance  peak  when  there  is  no  weak  signal  at  the  input  port  4  (A 
and  B).  When  a  weak  signal  is  input  at  port  4,  resonance  enhancement  is  observed  for  both 
the  unbroken-  and  broken-symmetry  cases  (C  and  D). 


When  the  gain  and  loss  are  balanced  but  the  coupling  strength  is  large  (distance  between  the 
resonators  is  small)  such  that  the  system  is  in  the  unbroken  PT-symmetric  phase  [Fig.  S9A(b) 
and  S9B(b)],  the  outputs  always  show  two  resonance  peaks  (mode  splitting)  regardless  of 
whether  the  field  is  input  at  port  1  or  port  3.  The  split  resonance  modes  have  the  same 
resonance  frequencies  and  similar  line  width  in  both  of  the  cases,  except  with  slight  variations 
in  their  heights. 
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When  the  gain  and  loss  are  balanced  and  the  coupling  strength  is  weak  (distance  between  the 
resonators  is  large)  such  that  the  system  is  in  the  broken  PT-symmetric  phase  [Fig.  S9A(c) 
and  S9B(c)],  only  the  output  port  2  shows  a  resonance  peak  regardless  of  whether  the  input  is 
at  port  1  or  at  port  3.  No  resonance  peak  or  dip  is  observed  in  the  output  port  4  in  both  cases. 
This  implies  that  the  input  field  is  localized  in  the  active  resonator  which  is  directly  coupled 
to  the  fiber  taper  with  the  ports  1  and  2. 


B8.  Effect  of  the  weak  probe  signal  on  the  output  spectra.  In  Fig.  S10,  we  give  enlarged 
spectra  of  those  shown  in  Figs.  4[b(ii)  and  b(iii)]  of  the  main  text.  It  is  clear  that  when  there 
is  no  weak  signal  at  the  input  port,  the  output  port  remains  at  the  noise  level  with  no 
observable  resonance  dip  or  peak  because  our  system  is  driven  below  the  lasing  threshold. 
One  can  say  that  the  photons  emitted  by  the  excited  ions  are  consumed  by  the  losses  (intrinsic 
and  coupling  losses)  in  the  system.  With  the  weak  signal  at  the  input  port,  we  see  a  resonance 
enhancement  at  the  output  port  for  both  the  broken-  and  the  unbroken-symmetry  cases. 


B9.  Supplementary  experimental  data  for  nonreciprocal  light  transmission.  For  the 
demonstration  of  nonreciprocal  light  transmission  in  our  PT-symmetric  system,  we  used  the 
configuration  given  in  Fig.  4S(B).  We  defined  the  transmission  from  port  1  to  port  4  as  the 
forward  T^4  and  transmission  from  port  4  to  port  1  as  the  backward  T4^i  transmission.  In 
the  main  text,  we  provided  the  results  for  the  T |  ,4  and  T4^,  demonstrating  that  in  the 
broken-symmetry  region,  forward  transmission  reduces  to  zero  T i _^4 ~  0  [Fig.  4a(iii)]  but  the 
backward  transmission  remains  high  [Fig.  4b(iii)].  In  the  experiments,  we  monitored  ports  2 
and  4  when  the  input  was  at  port  1 ,  and  ports  1  and  3  when  the  input  signal  was  at  port  4 .  In 
Fig.  4  of  the  main  text,  we  provided  only  the  spectra  at  port  4(1)  when  the  signal  input  was  at 
port  1  (4).  In  Fig.  Sll  below,  we  provide  the  transmission  spectra  at  ports  2  (3)  when  the 
input  signal  is  at  port  1  (4)  to  give  a  clear  understanding  of  the  signal  at  the  other  ports.  In  the 
broken-symmetry  region,  for  forward  propagation  we  observe  amplified  signal  at  output  2  at 
resonance  while  the  signal  at  output  4  is  almost  zero.  For  backward  propagation,  the 
amplified  signal  is  at  port  1;  however,  the  signal  at  port  3  has  a  shallow  resonance  dip. 


19 


NATURE  PHYSICS  |  www.nature.com/naturephysics 


©  2014  Macmillan  Publishers  Limited.  All  rights  reserved. 


SUPPLEMENTARY  INFORMATION 


DOI:  10.1038/NPHYS2927 


BIO.  Fluctuations  in  the  heights  of  resonance  peaks.  In  Fig.  3c  (unbroken  phase)  of  the 
main  text  and  in  Fig.  Sll,  the  heights  of  the  doublets  (peaks)  for  the  backward  and  forward 
transmission  show  a  difference.  This  difference  is  attributed  to  a  number  of  different  issues  in 
our  system.  First,  the  frequency  and  the  output  power  of  our  lasers  are  not  actively  controlled 
or  stabilized.  During  the  scanning  of  the  laser  frequency  to  obtain  the  transmission  spectra, 
the  outputs  of  the  lasers  fluctuate  which  is  reflected  as  the  fluctuations  at  the  heights  of  the 
peaks  or  dips  at  the  resonance  frequency.  Second,  since  we  do  not  have  active  stabilization 
and  control  mechanisms,  there  are  fluctuations  such  as  thermal  fluctuations  which  swing  the 
resonances  slightly  changing  the  measured  output  signal.  In  Fig.  S12,  we  provide  a  set  of 
experimental  results  obtained  on  the  same  setup  under  the  same  conditions  but  at  different 


Zero  gain 


Unbroken  symmetry  Broken  symmetry 


Fig.Sll.  Experimentally  obtained  spectra  in  broken-  and  unbroken-symmetry  regions.  (A) 

Output  at  port  2  when  the  input  signal  is  at  port  1  for  passive  photonic  molecule  [A(a)],  in 
PT-symmetric  system  when  the  symmetry  is  unbroken  [A(b)]and  when  it  is  broken  [A(c)]. 
Output  at  port  4  is  given  in  Fig.  4a  of  the  main  text.  (B)  Output  at  port  3  when  the  input 
signal  is  at  port  4  for  passive  photonic  molecule  [B(a)],  in  PT-symmetric  system  when  the 
symmetry  is  unbroken  [B(b)]and  when  it  is  broken  [B(c)] .  Output  at  port  1  is  given  in  Fig.  4b 
of  the  main  text. 
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times.  As  it  is  seen,  due  to  the  above  issues,  the  heights  of  the  peaks  fluctuate;  they 
sometimes  have  the  same  heights  but  at  some  other  times  one  is  higher  than  the  other. 

The  crucial  thing  to  point  out  here  is  that  even  if  there  are  fluctuations  in  the  heights  of  the 
peaks  for  backward  and  forward  transmissions  in  the  broken  and  unbroken  phases  in  the 
linear  regime  (Fig.  3c  of  the  main  text),  there  is  no  complete  absence  of  transmission  peak  in 
one  direction.  The  peaks  always  stay  there  regardless  of  the  fluctuations.  Note  that  when  the 
system  is  driven  to  the  nonlinear  regime,  for  the  broken  phase,  we  see  complete  absence  of 
transmission  peak  in  the  forward  direction  (Fig.  4a  of  the  main  text).  Similarly,  even  if  there 
are  slight  fluctuations  in  the  spectra  obtained  at  different  times  in  the  broken  phase  in  the 
nonlinear  regime  for  forward  and  backward  transmission,  there  is  always  a  strong 
transmission  in  one  direction  and  complete  absence  of  the  peak  in  the  other  direction.  Such 
fluctuation  can  be  minimized,  if  not  completely  eliminated,  by  active  stabilization  and  control 
of  laser  frequencies  and  output  powers  and  thermal  fluctuations  in  the  environment,  etc. 
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Fig.S12.  Transmission  spectra  obtained  at  different  times  when  the  system  is  in  the  linear 
regime.  Although  the  heights  of  the  resonance  peaks  fluctuate  in  measurements  taken  at 
different  times,  there  are  always  two  resonance  peaks  ( split  modes)  in  the  backward 
transmission  (red  spectra)  and  forward  transmission  (blue  spectra  in  the  inset)  directions 
when  the  system  is  in  the  linear  regime  and  PT -symmetry  is  unbroken. 
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Bll.  Comparison  of  nonreciprocity  demonstrated  in  PT -symmetric  microresonators  with 
reported  experiments  utilizing  nonlinearity  or  resonators.  Many  components  or  systems 
used  in  optics  are  reciprocal,  i.e.,  light  can  be  transmitted  in  both  directions.  Nonreciprocal 
devices  are  of  great  importance  and  much  needed  for  optical  communication  and  optical 
signal  processing.  For  example,  isolators  are  used  to  protect  the  laser  sources  and  sensitive 
components  from  back-reflected  light;  circulators  are  used  to  separate  and  route  light  in 
bidirectional  systems.  The  ability  to  control  the  direction  of  light  flow  in  such  a  way  that  light 
is  transmitted  in  one  direction  but  blocked  in  the  other  direction  requires  breaking  reciprocity 
or  the  time -reversal  symmetry.  In  many  optical  systems  used  today,  this  is  achieved  using 
magneto-optical  effects  induced  by  applied  magnetic  fields.  Unfortunately,  magneto-optic 
effect  in  many  materials  is  very  weak.  As  a  result,  magneto-optic  materials  with  large  sizes 
and  high  magnetic  fields  are  needed.  These  make  the  systems  complicated  and  bulky.  As  the 
technology  progresses  drastically,  there  is  a  drive  and  tendency  to  make  the  systems  smaller 
and  smaller,  and  if  possible  to  achieve  on-chip  optical  processing  systems  with  nano-  or 
micro-scale  footprints.  The  absence  of  magneto -optic  effects  in  materials  used  in 
conventional  optoelectronics  processing  demands  that  materials  with  higher  magneto-optic 
effect  are  suitably  integrated  into  the  on-chip  nano-or  micro-scale  structures.  However,  this  is 
not  an  easy  task.  Therefore,  there  is  an  ever-increasing  need  for  non-magneto-optic 
approaches  to  achieve  nonreciprocal  light  transmission.  Today  it  is  well-known  that 
reciprocity  can  be  broken  in  magneto-optic  materials,  nonlinear  materials  and  materials 
whose  dielectric  permittivity  and  magnetic  permeability  depend  on  time.  In  other  words,  a 
linear  static  dielectric  system  cannot  have  nonreciprocal  response11'12.  This  is  true  even  when 
gain  and  loss  exist  in  the  system1 1 .  Therefore,  a  system  without  magneto -optics  has  to  rely  on 
either  nonlinearity  or  time -dependent  effects  to  break  time -reversal  symmetry.  In  the 
following,  we  will  introduce  some  of  the  beautiful  experiments  demonstrating  nonreciprocal 
light  transmission  and  compare  their  performance  with  the  nonreciprocity  achieved  in  our 
work,  which  brings  PT-symmetric  concepts  with  nonlinearity  in  coupled  optical 
microresonators.  Our  focus  is  on  demonstrations  of  nonreciprocity  in  on-chip  waveguides 
and  microresonators,  and  nonreciprocity  based  on  nonlinear  effects  in  different  structures. 

Blla.  Nonreciprocity  based  on  nonlinear  effects.  Despite  many  theoretical  proposals,  there 
exist  only  a  very  small  number  of  experimental  demonstrations  of  nonreciprocal  light 
transmission  based  on  nonlinear  effects  such  as  second-order  nonlinearity  /2),  Kerr- 
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nonlinearity  x(3)  and  parametric  nonlinearity.  Raman  amplification  in  silicon  waveguides13'14, 
stimulated  Brillouin  scattering  in  silica  fiber15  and  dispersion-engineered  chalcogenide 
waveguides16,  %(2)  in  asymmetric  waveguides17  and  periodically  poled  lithium  niobate 
waveguide18,  and  a  nonlinear,  asymmetric,  distributed  Bragg  reflector  made  from  multilayer 
nonlinear  thin  film19,  stimulated  inter-polarization  scattering  in  a  photonic  crystal  fiber20,  and 
thermal  nonlinearity  in  a  pair  of  silicon  microring  resonators21,  just  to  name  a  few.  Among 
these  there  are  three  experimental  demonstrations  ’  ’  that  we  will  briefly  summarize 
below. 

Gallo  et  al.18  demonstrated  nonreciprocal  behavior  relying  on  the  tailoring  of  second  order 
nonlinearity  x(2)  with  quasi  phase  matching  gratings  in  a  LiNbC>3  channel  waveguide.  The 
gratings  were  4  cm  long  and  had  a  period  of  14.7  pm.  Moreover,  at  each  end  of  the  channel 
waveguide  1-mm-long  mode  filters  and  1.5-mm-long  mode  tapers  were  fabricated.  In  this 
device,  nonreciprocity  was  observed  for  input  peak  powers  beyond  1.5  W  at  A=  1 .55pm.  The 
contrast,  defined  as  the  difference  between  backward  and  forward  transmissions  normalized 
by  the  sum  of  the  transmissions  in  both  directions  at  the  wavelength  of  interest,  reached 
values  as  high  as  C=0.9  at  an  input  power  of  P=3.1  W.  Insertion  loss  of  the  device  was  ~0.7 
dB/cm. 

Kang  et  al.20  showed  nonreciprocal  light  transmission  in  a  photonic  crystal  fiber  (length  of 
1.5  m  and  a  core  of  1.8pm).  Here  breaking  the  reciprocity  was  achieved  via  the  stimulated 
inter-polarization  scattering  (SIPS)  which  is  a  nonlinear  process  relying  on  optical  excitation 
of  a  gigahertz  guided  acoustic  mode.  In  SIPS,  photons  from  a  pump  wave  travelling  in  the 
forward  direction  scatter  into  the  co-propagating  but  orthogonally  polarized  Stokes  wave  via 
the  excited  guided  acoustic  resonance.  Thus  in  the  forward  direction  the  pump  is  strongly 
attenuated.  However,  if  the  pump  is  launched  in  the  backward  direction  and  counter- 
propagates  with  the  Stokes  wave,  then  SIPS  does  not  take  place  and  the  pump  is  transmitted 
with  a  slight  attenuation  due  to  the  fiber  transmission  loss.  An  extinction  ratio  of  20  dB  was 
demonstrated  for  powers  of  ~  0.31  W. 

Fan  et  al.21 ,  on  the  other  hand,  used  a  system  of  two  silicon  microring  resonators  with  quality 
factors  of  Q~  27,000  and  Q~  43,800  configured  in  an  add-drop  filter  geometry  and 
asymmetrically  coupled  to  two  optical  waveguides  to  demonstrate  nonreciprocal  light 
transmission.  In  silicon  resonators,  two  photon  absorption  (TPA),  x(3)  and  thermal 
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nonlinearity  are  among  the  processes  that  can  lead  to  the  observed  nonreciprocal  behavior. 
The  authors  claimed  that  in  their  system  thermal  nonlinearity  was  the  main  reason  for  the 
observed  behavior.  As  a  result  with  input  powers  of  P=85  //W  and  P=850  p  W,  nonreciprocal 
transmission  ratios  of  ~20  dB  and  ~27  dB  were  demonstrated.  The  forward  insertion  loss  of 
the  device  was  12  dB. 

Bllb.  Nonreciprocity  in  optical  microresonators.  In  addition  to  the  experimental  work  by 
Fan  et  al.  where  nonreciprocity  relies  on  thermal  nonlineaity  in  silicon  microring  resonators 
and  our  present  work  which  relies  on  nonlinearity  enhancement  in  broken  PT-symmetry 
phase  for  nonreciprocal  light  transmission,  there  are  two  more  experimental  works  related 
with  optical  microresonators  worth  mentioning  here:  The  works  by  Bi  et  air  and  Tien  et 
al.23  where  magneto-optical  materials  are  integrated  with  microresonators.  We  should 
mention  that  there  are  many  proposals  on  how  nonreciprocity  can  be  achieved  by  integrating 
magneto-optical  materials  into  microresonators,  waveguides  and  photonic  crystal  structures 
either  using  sputter  deposition”  ,  wafer  bonding  or  adhesive  bonding”  .  The  concept  of 
nonreciprocal  ring  microresonator  was  introduced  by  Kono  et  al27  to  reduce  the  size  of 
isolators  and  other  nonreciprocal  devices  down  to  several  tens  of  micrometers  and  to  achieve 
on-chip  nonreciprocal  devices.  Here  we  will  introduce  only  the  works  by  Bi  et  air 2  and  Tien 
et  al23  which  rely  on  frequency  splitting  upon  the  application  of  magnetic  field  to  silicon 
microresonators  with  the  integrated  cerium-substituted  yttrium  iron  garnet  (Ce:YIG)  films. 
Whispering  gallery  mode  resonators  or  the  ring  resonators  support  two  counter-propagating 
modes  (clockwise  CW  and  counter-clockwise  CCW)  at  the  same  frequency.  An  magnetic 
field  applied  to  the  magneto -optical  microresonator  breaks  the  symmetry  of  the  resonator  and 
lifts  the  frequency  degeneracy22'23,28.  As  a  result,  the  CW  and  CCW  modes  have  different 
propagation  constants  and  resonance  frequencies,  thus  leading  to  non-reciprocal  optical 
transmission  in  the  wavelength  range  near  the  cavity  resonance. 

In  their  work23  Bi  et  al.  integrated  a  100  nm  garnet  film  layer  Ce:YIG(80  nm)/YIG(20  nm) 
using  a  pulsed  laser  deposition  technique  onto  a  silicon  racetrack  resonator.  They  then  lifted 
the  degeneracy  of  CW  and  CCW  resonant  modes  by  applying  a  magnetic  field  perpendicular 
to  the  direction  of  light  propagation.  An  extinction  of  19.5  dB  was  achieved  by  a  magnetic 
field  strength  of  1,500  Oe. 

Tien  et  al 23  integrated  Ce-YIG  of  thickness  500  nm  on  the  top  of  a  silicon  microring 
resonator  of  radius  900  pm  using  oxygen  plasma  enhanced  bonding.  They  achieved  an 
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extinction  of  9  dB  in  the  1550  nm  wavelength  band  with  an  applied  magnetic  field  strength  of 
50  Oe. 

Bile.  Other  interesting  methods  to  achieve  nonreciprocity.  Among  many  other  interesting 
approaches  to  achieve  nonreciprocal  light  transmission,  the  works  relying  on  opto- 
mechanical"  "  effects  and  interband  photonic  transitions  '  are  worth  mentioning  here. 

In  their  work  Manipatruni  et  al29  showed  that  optomechanical  systems  can  exhibit 
nonreciprocity  via  momentum  exchange  during  the  interaction  of  the  light  field  with  the 
mechanical  structure.  The  nonreciprocity  originates  from  the  asymmetry  of  the  radiation 
pressure  induced  by  the  photons  in  the  light  field  on  a  movable  mirror  for  incident  light  in 
opposite  directions.  They  then  propose  a  silicon  based  micro-optomechanical  device  that  can 
achieve  a  contrast  ratio  above  20  dB  for  an  incident  light  power  of  100  mW.  This  proposal  is 
yet  to  be  demonstrated  in  experiments. 

In  contrast  to  the  work  of  Manipatruni  et  al.29  where  nonreciprocity  is  based  on  a  nonlinear 
response  of  an  optomechanical  system,  Hafezi  and  Rabl30  proposed  another  optomechanical 
approach  which  can  lead  to  nonreciprocal  behavior  even  at  single  photon  level.  The  proposal 
relies  on  enhanced  nonlinear  coupling  between  a  mechanical  mode  of  a  WGM  microring 
resonator  and  an  optical  mode  with  a  probe  light  traveling  in  the  same  direction  as  the  strong 
pump  field.  This  leads  to  different  transmission  properties  for  the  probe  light  fields  traveling 
in  the  same  or  opposite  directions.  The  light  in  the  direction  of  the  pump  field  is  transmitted 
whereas  that  in  the  opposite  direction  is  blocked  by  the  resonator.  This  proposal  has  not  been 
realized  yet. 

In  2009,  Yu  and  Fan31  proposed  to  use  interband  transitions  induced  by  spatio-temporal 
refractive-index  modulations  in  photonic  structures  for  nonreciprocal  behavior  without 
magneto-optics.  The  nonreciprocal  behavior  relies  on  the  shifts  of  the  frequency  and 
wave  vector  during  the  photonic  transition.  In  a  structure  with  modulated  refractive  index,  the 
frequency  of  the  light  in  the  forward  direction  is  converted  to  a  higher  frequency  co2>co{\ 
however,  the  light  propagating  in  the  backward  direction  is  not  affected  by  the  modulation 
and  stays  intact.  If  an  absorption  filter  that  absorbs  the  light  at  co2  completely  is  placed  at  the 

output,  the  system  will  behave  as  an  isolator  which  blocks  the  light  in  the  forward  direction 
but  transmits  it  in  the  backward  direction.  An  interesting  property  of  this  approach  is  that  it  is 
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a  linear  process  and  that  the  photonic  transition  does  not  depend  on  the  phase  or  the 
amplitude  of  the  incident  light. 

Although  Yu  and  Fan11  proposed  to  use  a  microring  resonator  for  experimental  realization,  in 
2012  Lira  et  al.32  demonstrated  their  proposal  and  achieved  nonreciprocity  in  a  silicon  slotted 
waveguide  structure  that  was  electrically  driven  to  induce  indirect  interband  photonic 
transition.  Here,  an  indirect  interband  photonic  transition  took  place  between  two  optical 
modes  having  different  longitudinal  wave  vectors,  and  transverse  modal  profiles  with 
different  symmetries.  By  applying  an  electrical  modulation  of  10  GHz,  corresponding  to  the 
frequency  difference  between  the  optical  modes  of  interest  in  their  slotted  waveguide,  the 
authors  achieved  a  contrast  of  3  dB  at  an  applied  electrical  signal  power  of  25  dBm.  The 
insertion  loss  of  the  system  was  around  70  dB . 

Blld.  What  makes  our  scheme  different  than  above  mentioned  experimental  realizations 
and  what  are  the  advantages?  In  our  work,  we  used  two  directly  coupled  microtoroid 
resonators  (quality  factors  of  the  order  of  107)  configured  in  an  add-drop  filter  structure.  One 
of  the  resonators  is  an  active  resonator  with  erbium  ions  as  the  embedded  gain  dopants  within 
silica  matrix  whereas  the  other  has  passive  silica  loss  at  the  k=l  ,55pm  band.  Although  the 
demonstrated  nonreciprocity  in  our  work  also  relies  on  nonlinearity,  conceptually  it  is  very 
different  than  the  above  works.  The  key  point  in  our  work  is  the  use  of  PT-symmetric 
concept,  where  the  loss  in  passive  resonator  is  balanced  with  the  gain  in  the  active  resonator 
and  the  coupling  strength  between  the  resonators  is  adjusted  such  that  the  system  operates  in 
the  broken  PT  phase.  As  a  result,  in  the  broken  PT-symmetry  phase  the  optical  field  is 
strongly  localized  in  the  resonator  with  gain,  which  in  turn  enhances  the  nonlinear  process 
(i.e.,  the  nonlinearity  can  be  observed  with  low  power  levels).  In  the  unbroken  phase  there  is 
no  field  localization  and  hence  one  needs  higher  power  levels  to  observe  the  nonlinearity. 
This  is  seen  clearly  in  Fig.  S9  and  Fig.  3a  of  the  main  text.  There  is  a  significant  difference 
between  the  observed  nonreciprocity  in  the  broken-  and  unbroken-phases. 

Up  to  date  there  was  no  reported  experiment  in  the  literature  which  utilizes  PT-symmetric 
structures  to  achieve  nonreciprocal  light  transmission.  Our  work  is  the  first  experimental 
demonstration  bringing  together  the  PT -symmetric  concepts  with  nonlinearity  to  demonstrate 
nonreciprocal  light  transmission  which  is  crucial  for  many  optical  devices  and  photonic 
applications.  The  advantages  of  our  scheme,  which  brings  together  PT-symmetric  concepts 
with  nonlinearity-induced  nonreciprocal  light  transmission,  over  the  non-PT  schemes 
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utilizing  nonlinearity  are  significant  reduction  in  the  input  power  to  observe  nonreciprocity 
(~lpW  in  this  work  versus  3W  in  Ref.  [18],  0.310W  in  Ref.  [20]  and  ~  85pW  in  Ref.  [21]), 
higher  contrast,  smaller  footprint  and  complete  absence  of  the  signal  in  one  direction  but 
resonantly  enhanced  transmission  in  the  other  direction.  Unlike  all  the  experimental  works 
mentioned  above  (those  utilizing  magneto-optical  effects,  nonlinearity  and  interband 
transitions),  we  observed  in  this  work  a  complete  absence  of  resonance  peak  in  one  direction. 
Our  work  constitutes  the  first  direct  experimental  proof  of  the  connection  between 
nonreciprocity  and  PT-symmetry  which  has  been  largely  confused  [see  the  work  by  Feng  et 
al ,33  and  comments  on  this  article  by  Fan  et  al.n,  as  well  as  the  paper  Wang  et  al?4  and  the 
comment  by  Petrov  et  al.35]. 


C.  Imperfect  gain/loss  balance  in  PT-symmetric  systems 

In  PT-symmetric  systems,  formed  by  coupled  structures  such  as  waveguides  or  resonators, 
with  exactly  balanced  gain  and  loss,  the  phase  transition  ( PT-symmetry  breaking)  from 
broken  to  unbroken  symmetry  is  a  sharp  bifurcation  in  orthogonal  directions  in  the  real  parts 
of  the  eigenfrequencies  of  the  supermodes.  This  is  accompanied  by  a  sharp  orthogonal 
coalescence  of  imaginary  parts  of  the  eigenfrequencies.  However,  in  our  experiments,  we 
observed  that  the  bifurcations  were  not  so  sharp  and  orthogonal  as  predicted  by  the  theory  of 
perfectly  balanced  systems,  but  rather  there  is  a  smooth  gradual  separation.  We  attributed  this 
to  the  imperfection  in  gain/loss  balance. 

To  explain  this  smooth  bifurcation,  we  formulate  a  theoretical  model  in  which  the  gain  and 
loss  are  not  perfectly  balanced36.  We  construct  the  equations  of  motion  of  linearly  coupled 
oscillators  x  and  y  , 


d2x  dx  2 

— -  +  u - 1-  co~x  +  k  y  =  0 

dt 2  dt 


dy  2  „ 

-v - \-co  y  +  icx  =  0 

dt 


(S.9) 
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Both  oscillators  have  the  same  natural  frequency  co ,  the  parameters  //  and  v  are  a  measure 
of  the  loss  and  the  gain,  and  k  is  the  coupling  strength.  We  seek  solutions  of  the  form  em , 
which  lead  to  the  quartic  polynomial  equation 

A4 -i(p-v)X3  -(2 co2  -pv^A2  +ico2  (jU-v)A-k1  +co4  =0  (S.10) 


Fig.S13.  Real  and  imaginary  parts  of  the  eigenfrequencies  of  the  coupled  system  as  a 
function  of  the  coupling  strength.  The  case  where  both  of  the  oscillators  are  lossy 
f  t  =  -v  =  0.01  is  depicted  in  green  diamond  marks  [A(a)  and  B(a) ].  Exactly  balanced  gain 
and  loss  cases  for  small  loss  ( p-v  =0.01 )  and  larger  loss  ( p  -v  -  0.02 )  are  shown  in  [A(b) 
and  B(b)]and  [A(c)  and  B(c)],  respectively,  with  brown  crosses  and  blue  solid  circles.  Real 
and  imaginary  parts  of  the  eigenfrequencies  for  unbalanced  loss  and  gain  are  given  in  (C) 
and  (D),  respectively.  Loss  and  gain  parameters  are  set  as  ( p  -  3v  =  0.03)  and  ( //  =  8u  =0.08 
),  respectively,  for  [C(a)  and  D(a)]  brown  cross  marks  and  for  [C(b)  and  D(b)]  blue  solid 
circle  marks. 
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We  have  numerically  solved  this  equation  for  co  =  1  at  various  coupling  strengths  and  gain-to- 
loss  ratios.  The  real  and  imaginary  parts  of  the  eigenfrequencies  of  this  coupled  system  are 
plotted  in  Fig.  S13  as  functions  of  coupling  strength  k  and  for  chosen  values  of  gain-to-loss 
ratios.  We  see  that  for  the  case  of  exactly  balanced  loss  and  gain,  the  bifurcation  of  the  real 
and  imaginary  parts  of  the  eigenfrequencies  at  the  phase  transition  point  is  sharp  and  in 
orthogonal  directions  (Figs.  S13A  and  S13B).  However,  for  the  unbalanced  case  the 
bifurcations  are  not  abrupt,  but  rather  are  smooth  (Figs.  S13C  and  S13D).  The  degree  of 
smoothness  increases  with  increasing  imbalance  between  gain  and  loss.  Moreover,  for  the 
unbalanced  loss-and-gain  case  the  eigenfrequencies  are  never  real  and  there  is  always  a 
nonzero  imaginary  part. 

It  is  also  of  critical  importance  to  notice  that  for  the  special  case  in  which  the  loss  and  gain 
are  exactly  balanced,  the  coupled  equations  of  motion  for  and  y(7)  above  can  be 
derived  from  a  Hamiltonian  H  , 


H  =  pq  +  y[yq-xp}  +  {a>2  -^2)xy  +  —  k(x1  +  y2) 


(S.ll) 


where  p  and  q  are  momenta  conjugate  to  x  and  y  ,  and  2 y  =  ju  =  v  .  (The  existence  of  a 
Hamiltonian  is  surprising  because  the  system  has  loss  and  gain.)  In  this  case,  the  energy  H 
of  the  system  is  exactly  conserved,  although  it  is  not  a  simple  expression  such  as  the  sum  of 
the  squares  of  the  momenta  and  the  coordinates.  Interestingly,  if  the  coupling  strength  k 
becomes  strong  enough,  the  frequencies  become  complex.  Thus,  there  are  two  regions  of 
broken  PT  symmetry,  one  for  weak  coupling  and  one  for  strong  coupling’6.  Because  the 
system  is  Hamiltonian,  it  can  be  quantized  by  imposing  the  requirement  that  x  and  p  (and 
also  y  and  q )  obey  equal-time  commutation  relations.  One  can  then  find  the  quantized 
energies  of  the  Hamiltonian.  One  obtains  the  remarkable  result  that  the  quantum  energies  are 
real  for  exactly  the  same  range  of  parameters  that  the  classical  frequencies  are  real;  that  is  the 
region  of  unbroken  PT  symmetry.  The  quantum  energies  become  complex  when  the  classical 
frequencies  are  complex;  that  is  in  the  region  of  broken  PT  symmetry36. 

Our  experimental  results  depicted  in  Fig.  2  of  the  main  text  show  the  same  behavior  as  the 
numerical  results  for  our  theoretical  model  shown  in  Fig.  S13  C(a)  and  D(a),  implying  that 
the  gain  and  loss  in  the  experiments  are  not  exactly  balanced.  In  Figs.  S13  A(a)  and  B(a),  we 
see  that  for  two  coupled  lossy  oscillators  (passive  resonators),  the  difference  of  the  real  parts 
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of  the  eigenfrequencies  increases  with  increasing  coupling  strength,  whereas  the  imaginary 
parts  remain  equal.  These  numerical  predictions  agree  well  with  the  results  of  our 
experiments  as  shown  in  Figs.  2a  and  b  of  the  main  text;  that  is,  the  mode  splitting 
(difference  in  the  real  parts  of  the  eigenfrequencies)  increases  with  increasing  coupling 
strength  while  the  difference  in  the  imaginary  parts  of  the  eigenfrequencies  stays  the  same. 
Finally,  as  depicted  in  Fig.  2c  and  d,  we  have  observed  in  the  experiments  that  the  higher  the 
initial  loss  (lower  Q)  of  the  resonators,  the  higher  the  coupling  strength  needed  to  observe 
the  transition  from  broken-symmetry  to  unbroken-symmetry  region.  This  is  indeed  what  is 
found  in  the  numerical  solutions  of  the  characteristic  equation  of  the  theoretical  model  (Fig. 
SI 3).  Thus,  the  theoretical  model  introduced  here  and  our  experimental  observations  are 
consistent.  Our  model  is  physically  realistic  because  in  practical  realizations  it  is  almost 
impossible  to  have  exactly  balanced  loss  and  gain. 

The  observation  in  our  experiments  and  the  theoretical  model  above  is  critical  in  two  ways. 
First,  imposing  exactly  balanced  loss  and  gain  in  coupled  systems  for  extended  durations  of 
time  is  not  practical.  Second,  In  the  broken-symmetry  phase,  the  field  propagating  in  the  PT- 
symmetric  system  is  always  confined  in  the  structure  with  gain,  thereby  experiencing  a  strong 
overall  gain  and  leading  to  enhanced  transmission.  At  the  phase  transition  point,  the  change 
in  the  gain  experienced  by  the  field  should  be  abrupt  in  the  ideally  balanced  gain  and  loss 
situation.  Thus,  if  the  coupling  strength  is  changed  by  a  very  small  amount  around  the  phase 
transition  point,  there  is  a  large  abrupt  change  in  the  real  and  imaginary  parts  of  the 
eigenfrequencies,  leading  to  an  abrupt  change  in  the  amplification  or  the  gain  experienced  by 
the  field.  In  the  non-ideal  case  where  the  bifurcation  is  not  orthogonal  but  smoothened,  the 
change  in  the  gain/amplification  experienced  by  the  field  around  the  phase  transition  point  is 
reduced. 

Finally,  we  should  note  here  that  the  smooth  bifurcation  issue  was  also  discussed  by  Benisty 
et  al ,37  from  a  different  perspective:  They  attributed  the  smoothness  around  the  phase 
transition  point  to  a  complex  coupling  arising  as  the  gain  is  tuned  in  one  of  the  systems. 
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